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We study the quasiparticles of the one-band t-J and t-t'-t"-J models using a variational ap¬ 
proximation that includes spin fluctuations in the vicinity of the hole. We explain why the spin 
fluctuations and the longer range hopping have complementary contributions to the quasiparticle 
dynamics, and thus why both are essential to obtain a dispersion in agreement with that measured 
experimentally. This is very different from the three-band Emery model in the strongly-correlated 
limit, where the same variational approximation shows that spin fluctuations have a minor effect 
on the quasiparticle dynamics. This difference proves that these one-band and three-band models 
describe qualitatively different quasiparticles, and therefore they they cannot both be suitable to 
describe the physics of underdoped cuprates. 
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I. INTRODUCTION 

Nearly three decades after their discovery [I|, the high- 
temperature cuprate superconductors have so far eluded 
a comprehensive explanation. These layered materials 
contain Cu02 layers, which are antiferromagnetic insu¬ 
lators in the undoped limit and become superconduct¬ 
ing upon doping Q. The hole-doped side shows a more 
robust superconductivity, extending to higher tempera¬ 
tures and over a wider range of dopings. The first step 
towards deciphering the mechanism of superconductivity 
in these compounds is a proper description of the mo¬ 
tion of these holes in the Cu02 layer - this has become 
one of the most studied problems in condensed matter 
theory ii- Its solution should elucidate the nature of 
the quasiparticles that eventually bind together into the 
pairs that facilitate the superconducting state. 

Despite significant effort, even what is the minimal 
model that correctly describes this low-energy quasipar¬ 
ticle is still not clear. There is general agreement that 
the parent compounds are charge-transfer insulators 
and wide consensus that most of their low-energy physics 
is revealed by studies of a single Cu02 layer, modeled in 
terms of Cu and O 2 p orbitals. Because only 

ligand 2 p orbitals hybridize with the ^d„.2_y2 orbitals, it 
is customary to ignore the other O 2 p orbital^ this leads 
to the well-known three-band Emery model (3. 

However, the Emery model is perceived as too compli¬ 
cated to study so it is often further simplified to a one- 
band t-J model that describes the dynamics of a Zhang- 
Rice singlet (ZRS) @,0. We known that the t-J model 
with only nearest-neighbor (nn) hopping t is certainly not 
the correct model because it predicts a nearly flat quasi¬ 
particle energy along (0 , tt) — (tt, Ob unlike the substantial 
dispersion found experimentally |8l-ll0l|. However, its ex¬ 
tension with longer range hopping, the t-t'-t"-J model, 
was shown to reproduce the correct dispersion m for 
values of the 2"'^ and 3'^'^ nn hoppings in rough agreement 
with those estimated from density functional theory [l^ 
and cluster calculations [l^. Taking this agreement as 


proof that this is the correct model, both the t-J model 
and its parent, the Hubbard model [3, have been stud¬ 
ied very extensively in the context of cuprate physics. 

Here we show that the quasiparticle of the t-t'-t''-J 
model is qualitatively different from that of the Uud 00 
limit of the three-band Emery model [l^. While both 
models predict a dispersion in quantitative agreement 
with that measured experimentally (for suitable values of 
the parameters), the factors controlling the quasiparticle 
dynamics are very different. It has long been known that 
both the longer-range hopping and the spin fluctuations 
play a key role in the dynamics of the quasiparticle of the 
t-t'-t"-J model. Here we use a non-perturbative varia¬ 
tional method, which agrees well with available exact di- 
agonalization (ED) results, to show that spin-fluctuations 
and longer-range hopping control the quasiparticle dis¬ 
persion in different parts of the Brillouin zone, and to 
explain why. In contrast, using the same variational 
approach, it was recently argued that spin fluctuations 
play no role in the dispersion of the quasiparticle of the 
Udd —>■ 00 limit of the Emery model [1^. This claim is 
supported by additional results we present here. 

This major difference in the role played by spin fluctu¬ 
ations in determining the quasiparticle dynamics shows 
that these models do not describe the same physics. This 
suggests that the t-t'-t''-J model is not suitable for the 
study of cuprates in the hole-doped regime, although it 
and related one-band models may be valid in the electron- 
doped regime. As we argue below, it may be possible to 
“fix” one-band models by addition of other terms, al¬ 
though we do not expect this to be a fruitful enterprise. 
Instead, we believe that what is needed is a concerted 
effort to understand the predictions of the Emery model. 
Our results in Ref. [T^ and here that spin fluctuations of 
the AFM background do not play a key role in the quasi¬ 
particle dynamics of this three-band model, contrary to 
what was believed to be the case based on results from 
one-band models, should simplify this task. 

The article is organized as follows. In Section H we 
review the three-band Emery model and briefly discuss 
the emergence of the one-band and simplified three-band 
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models in the asymptotic limit of strong correlations 
on the Cu sites. Section III describes the variational 
method, which consists in keeping a limited number of 
allowed magnon configurations in the quasiparticle cloud. 
Section IV presents our results for both one- and three- 
band models, and their interpretation. Finally, Section 
V contains a summary and a detailed discussion of the 
implications of these results. 

II. MODELS 

A widely accepted starting point for the description of 
a Cu02 layer is the three-band Emery model 

= Tpp “t” T^pd “b ^pd ^ ^ 

2^0,cr 

“t“ Upp E + Udd X! (1) 

The sets O and Cu contain the ligand O 2p and the Cu 
orbitals respectively, sketched in Fig. [TJa). For 
J € O, Tii^cr = pI is the number of spin-cr holes in that 
2p orbital. Similar notation is used for the 3d orbitals, 
their hole creation operators being d ■ ^, i G Cu. 

Tpp is the kinetic energy of the holes moving on the 
O sublattice, described by a Hamiltonian with first (tpp) 
and second (tpp) nearest-neighbour (nn) hopping; 

Tpp = tpp 'y ^ ^Sp\^cTPi+d,(T 

iGO,S,a- 

^PP ^ ^ Pi^aiPi—e.crTPiJ^^ u)- (2) 

The lattice constant is set to a = 1. The vectors 6 = 
±(0.5, 0.5), ±(0.5,—0.5) are the distances between any 
O and its four nn O sites, and = ±1 sets the sign 
of each nn pp hopping integral in accordance with the 
overlap of the 2p orbitals involved, see Fig. [TJa). Next nn 
hopping is included only between O 2p orbitals pointing 
toward a common bridging Cu, separated by e = (1,0) 
or (0,1); hybridization with the 4s orbital of the bridging 
Cu further boosts the value of this hopping integral. 

Tpd is the kinetic energy of holes moving between 
neighbour Cu and O orbitals: 

Tpd = tpii y ( ^ud| ± h.c., (3) 

2GCu,U,(7 

where u = (±0.5, 0), (0, ±0.5) are the distances between 
a Cu and its four nn O sites, and are the signs of 
the overlaps of the corresponding orbitals. It is this term 
that provides the main justification for ignoring the other 
sets of 2p orbitals, because symmetry forbids hopping 
of Cu holes from orbitals into the non-ligand O 

orbitals. We further discuss this assumption below. 

Apt; is the charge transfer energy which ensures that in 
the parent compound the O 2p orbitals are fully occupied 


(a) (b) 



FIG. 1. (color online) (a) Sketch of the Cu02 layer, with O are 
marked by squares and Cu marked by circles. The relevant 
orbitals are drawn at a few sites. The arrows indicate the 
hopping terms included in the three-band Emery model; (b) 
Sketch of the one-band t-t'-t"-J model. Cu sites host spin 
degrees of freedom, except at sites where a Zhang-Rice singlet 
is centered (red circle). The arrows indicate the various terms 
in this Hamiltonian; (c) Sketch of the Udd —>■ oo limit of the 
Emery model. Cu sites host spin degrees of freedom but the 
doped holes (red filled square) move on the O lattice. The 
arrows indicate the various terms in this Hamiltonian; (d) 
The full Brillouin zone of the Cu02 lattice (the outer square) 
which encloses the magnetic Brillouin zone (shaded). 

(i.e. contain no holes). Finally, Upp and Udd are the 
Hubbard repulsion in the 2p and 3d orbitals, respectively. 
Longer range Coulomb interaction between holes on O 
and Cu can also be added, but for the single doped-hole 
problem analyzed here, it leads to a trivial energy shift. 

Strong correlations due to the large Udd, combined 
with the big Hilbert space with its three-orbitals basis, 
and the need for a solution for a hole concentration equal 
(in the parent compound) or larger (in the hole doped 
case) than one per Cu, make this problem very difficult 
to solve. While progress is been made with a variety of 
techniques [13 (which however have various restrictions, 
such as rather high-temperatures and/or small clusters 
for quantum Monte Carlo methods, and/or additional 
approximations, such as setting tpp = 0 for convenience), 
it is far more customary to further simplify this model 
before attempting a solution. 

A reasonable way forward is to use the limit Udd oo 
to forbid double occupancy of the 3d orbitals. As a re¬ 
sult, in the undoped ground-state there is one hole per 
Cu. Virtual hopping processes lead to antiferromagnetic 
(AFM) superexchange between the resulting spin degrees 
of freedom, so that the parent compound is a Mott insu¬ 
lator with long-range AFM order |la |. 

Upon doping, holes enter the O band and the issue is 
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how to accurately describe their dynamics as they inter¬ 
act with the spins at the Cu sites. Here we compare two 
such descriptions for the single doped hole case. 


A. One-band models 


In their seminal work Q , Zhang and Rice argued that 
the doped hole occupies the linear combination of the 
four O 2p ligand orbitals surrounding a central Cu, that 
has the same symmetry like the Cu 3d orbital 

hosting the spin. Furthermore, exchange locks the two 
holes in a low-energy Zhang-Rice singlet (ZRS). They 
also argued that the dynamics of this composite object, 
which combines charge and spin degrees of freedom, is 
well captured by the one-band model: 


n = r 


T + T' + T‘ 


V -t- T^afm- 


( 4 ) 


The first term describes the hopping of the ZRS (marked 
in Fig. mb) by the “missing spin” locked in the ZRS) 
on the square lattice of Cu sites that hosts it. Originally 
only nn hopping T was included: T = + 

/i.c., with i,j € Cu. The projector V removes doubly 
occupied states, therefore this term allows only Cu spins 
neighboring the ZRS to exchange their location with the 
ZRS. This mimics the more complex reality of the doped 
hole moving on the O sublattice and forming ZRS with 
different Cu spins. 

Although in Ref. @ it was argued that only nn ZRS 
hoping is important, longer-range 2nd (T') and 3rd (T") 
nn hopping was later added the model on a rather ad-hoc 
basis. As discussed below, this is needed in order to find 
a quasiparticle dispersion similar to that measured exper¬ 
imentally. These terms are defined similarly to T with 
hopping integrals t' and t", respectively. For cuprates, 
t'/t ^ —0.3, t”/t ~ 0.2 are considered to be representa¬ 
tive values , in agreement with various estimates 

[I10. In the following, we refer to this as the 
model, whereas if t' = t" = 0 we call it the t-J model. 

The term 'Hafm = describes the nn 

AFM superexchange between the other Cu spins Si, with 
J/t ~ 0.3 for cuprates [2^. It leads to AFM order in the 
undoped system [T^, and also controls the energy of the 
cloud of magnons that are created in the vicinity of the 
ZRS, as it moves through the magnetic background. 

The t-J model also emerges as the C/ —>■ oo limit of the 
Hubbard model 0 , but with additional terms of order 
J. One of them, —J/4,J2{ij)''^inj, gives trivial energy 
shifts for both the undoped and the single-hole doped 
cases of interest to us in this work, so its presence can 
be safely ignored in this context. More interesting is the 
so-called three-site term VT^s'P [HI, where 
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describes ZRS hopping through an intermediate Cu site 
and permits spin swapping with the spin at this inter¬ 
mediate site. As shown below, this term influences the 
quasiparticle dispersion but it is not clear that it should 
be included in the one-band model, because the original 
Hamiltonian is the Emery, not the Hubbard, model. 

In fact, a perturbational derivation of the low-energy 
Hamiltonian obtained by projecting the three-band 
model onto ZRS states reveals a much more complicated 
Hamiltonian than the t-t'-t"-J model, with many other 
terms 0. We are not aware of a systematic study of 
their impacts, but their presence underlies one important 
issue with this approach: the hoped-for simplification due 
to the significant decrease in the size of the Hilbert space 
comes at the expense of a Hamiltonian whose full expres¬ 
sion is very complicated. Using instead simpler ver¬ 
sions like the t-t'-t''-J model may result in qualitatively 
different physics than that of the full one-band model. 
Here we argue that this is indeed the case. 


B. Simplifled three-band model 

An alternative is to begin at the same starting point, 
i.e. the limit Udd oo resulting in spin degrees of free¬ 
dom at the Cu sites. However, the O sublattice on which 
the doped hole moves is kept in the model, not projected 
out like in the one-band approach, see Fig. [T])c). This 
leads to a bigger Hilbert space than for one-band models 
(yet smaller than for the Emery model) but because spin 
and charge degrees of freedom are no longer lumped to¬ 
gether, the resulting low-energy Hamiltonian is simpler 
and makes it easier to understand its physics. 

The effective model for a layer with a single doped hole, 
which for convenience we continue to call ’’the three-band 
model” although it is its Udd ^ oo approximation, was 
derived in Ref. 0 and reads: 

Hes = 'Hafm + -\- ^pp ^swap ■ (6) 

The meaning of its terms is as follows: 

'Hafm = Jdd ^ • Sj 

(iJ)' 

is again the nn AFM superexchange between the Cu spins 
Si, so Jdd = J oi the one-band models. The main dif¬ 
ference is indicated by the presence of the prime, which 
reflects the absence of coupling for the pair that has the 
doping hole on their bridging O. The next term, 

T~^Jpd ~ ’^pd 'y ^ Sj ■ Sj-fu 
ieO,u 


is the exchange of the hole’s spin = 4 '^a,jdP\a^apPip 
with its two nn Cu spins. It arises from virtual hopping 
of a hole between a Cu and the O hosting the doped hole. 

Like in Eq. m.Tpp is the kinetic energy of the doping 
hole as it moves on the O sublattice. It is supplemented 
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by Tswap which describes effective hopping mediated by 
virtual processes where a Cu hole hops onto an empty 
O orbital, followed by the doping hole filling the now 
empty Cu state [2^. This results in effective nn or next 
nn hopping of the doped hole, with a swapping of its spin 
with that of the Cu involved in the process. The explicit 
form of this term is: 


T — —f 

swap — ^ 


-EE 

iGCu,U:^u' ( 7 ,( 7 ' 


Su-1 


Pi-\-u,aPi-\-u' ,(7' 




reflecting the change of the Cu spin located at from 
a to a' as the hole changes its spin from a' to a while 
moving to another O. The sign = ±1 is due to the 
overlaps of the 2p and 3d orbitals involved in the process. 

For typical values of the parameters of the Emery 
model 0 and using Jm (~ 0.15 eV) as the unit of en¬ 
ergy, the dimensionless values of the other parameters 
are tpp ~ 4.1, t'pp ^ OMpp, ^ 3.0 and Jpd ^ 2.8. We 
use these values in the following, noting that the results 
are not qualitatively changed if they are varied within 
reasonable ranges. For complete technical details of the 
derivation of this effective Hamiltonian and further dis¬ 
cussions of higher order terms, as well as a comparison 
with other work along similar lines [2^, the reader is 
referred to the supplemental material of Ref. [2^ . 


III. METHOD 

The ground state of the undoped layer is not a simple 
Neel-ordered state. This is due to the spin fluctuations 
term = J/2 i^r E + S ~) present in "Hafm, 

which play an important role in lower dimensions. A 2D 
solution can only be obtained numerically, for finite size 
systems [2^. The absence of an analytic description of 
the AFM background has been an important barrier to 
understanding what happens upon doping, because the 
undoped state itself is so complex. It is also the reason 
why most progress has been computational in nature and 
mostly restricted to finite clusters. While such results are 
very valuable, it can be rather hard to gauge the finite- 
size effects and, more importantly, to gain intuition about 
the meaning of the results. 

Because our goal is to verify whether the two kinds 
of models have equivalent quasiparticles, which requires 
us to understand qualitatively what controls their dy¬ 
namics, we take a different approach. We use a quasi- 
analytic variational method valid for an infinite layer, so 
that finite-size effects are irrelevant. By systematically 
increasing the variational space we can gauge the accu¬ 
racy of our guesses and, moreover, also gain intuition 
about the importance of various configurations and the 
role played by various terms in the Hamiltonians. Where 
possible, we compare our results with those obtained by 
exact diagonalization (ED) for small clusters, providing 
further proof for the validity of our method. 

Eor simplicity, in the following we focus on the one- 
band model; the three-band model is treated similarly, as 


already discussed in Ref. [I^ . Because we do not have an 
analytic description of the AEM background wavefunc- 
tion, we divide the task into two steps. 


A. Quasiparticle in a Neel background 


In the first step we completely ignore the spin fluctua¬ 
tions by setting 'Hsf 0, to obtain the so-called 
model. As a result, the undoped layer is described by a 
Neel state |N) with up/down spins on the A/B sublat¬ 
tice, without any spin-fluctuations. One may expect this 
to be a very bad starting point, given the importance of 
spin-fluctuations for a 2D AFM. At the very least, this 
will allow us to gauge how important these spin fluctua¬ 
tions really are, insofar as the quasiparticle dynamics is 
concerned, when we include them in step two. 

It is also worth remembering that the cuprates are 
3D systems with long-range AFM order stabilized up to 
rather high temperatures by inter-layer coupling, in the 
undoped compounds. The spin fluctuations must there¬ 
fore be much less significant in the undoped state than 
is the case for a 2D layer, so our starting point may be 
closer to reality than a wavefunction containing the full 
description of the 2D spin fluctuations. 

Magnons (spins wrongly oriented with respect to their 
sublattice) are created or removed when the ZRS hops 
between the two magnetic sublattices. The creation of 
an additional magnon costs up to 2J in Ising exchange 
energy as up to four bonds involving the magnon now 
become FM. This naturally suggests the introduction of 
a variational space in terms of the maximum number of 
magnons included in the calculation. 

This variational calculation is a direct generalization 
of that of Ref. [2^, where the quasiparticle of the t-J^ 
model was studied including configurations with up to 7 
adjacent magnons. That work showed that keeping con¬ 
figurations with up to three magnons is already accurate 
if t/ J is not too large, so here we restrict ourselves to this 
smaller variational space. (Note that three is the mini¬ 
mum number of magnons to allow for Trugman loops 27] , 
see discussion below, so a lower cutoff is not acceptable). 
The configurations included are the same as in Ref. , 
where this type of approach was first pioneered. 

To be more specific, our goal is to calculate the one- 
hole retarded Green’s function at zero temperature: 


G(k,a;) = (N|4,G(cc)du.t|N), 


( 7 ) 


where G(w) = lim,j_>Q+(a; — % + irf) ^ is the re¬ 
solvent of the one-band Hamiltonian o and dk.t = 
Here A —> oo is the number of 
sites in each magnetic sublattice, k is restricted to the 
magnetic Brillouin zone depicted in Fig. [IJd), and we 
set = I. The spectrum is identical if the quasiparticle 
is located on the spin-down sublattice: conservation of 
the total z-axis spin guarantees that there is no mixing 
between these two subspaces with different spin. 
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Taking the appropriate matrix element of the identity 
G{oj){uj + ir] — H) = 1 leads to the equation of motion; 

[uj + ir] — J — e(k)]G(k, w) — t -P'i(k, w, e) = 1. (8) 

e 

The four vectors e = ±(1, 0), ±(0,1) point to the nn Cu 
sites, J is the Ising exchange energy cost for adding the 
ZRS (four AFM bonds are removed) and 

e(k) = At' cos kx cosky + 2t"[cos{2kx) + cos(2/cy)] (9) 

is the kinetic energy of the ZRS moving freely on its 
own magnetic sublattice. NN hopping creates a magnon 
as the hole moves to the other magnetic sublattice; this 
introduces the one-magnon propagators: 

Ai(k,a;,e) = ^ ^ e*‘^'^‘(N|4^G(c.)5-d.+,.jN) 

ieCuA 

with the hole on the B sublattice and therefore a present 
magnon on a nn A site, to conserve the total spin. 

Equation ([5]) is exact (for an Ising background) but 
it requires the Fi(k,a;,e) = Fi{e) propagators to solve. 
Their equations of motion (EOM) are obtained similarly: 

[cu + iT^- ^J]Fi(e) = t'J2 

e'-Le 

+ t[G(k, w) + E 2 (e, e) + E 2 (e, e')]. (10) 

e'_Le 

Note that 2'^'^ and 3“^ nn hopping keeps the hole on the 
B sublattice and thus conserve the number of magnons, 
linking Fi to other Fi propagators. However, nn hop¬ 
ping links Fi to G(k,a;) if the hole hops back to the 
A sublattice by removing the existing magnon, but also 
to two-magnon propagators, F 2 , if it hops to a differ¬ 
ent A site than that hosting the first magnon. The 
equation above imposes the variational restriction that 
the two magnons are adjacent, so only F 2 (k, w,e, e') = 

S^GCu .4 with C-k 

e' ^ 0 are kept. This is a good approximation for 
the low-energy quasiparticle whose magnons are bound 
in its cloud, and thus spatially close. (Because fewer 
AFM bonds are disrupted, these configurations cost less 
exchange energy than those with the magnons apart). 
Of course, the hole could also travel far from the first 
magnon (using 2'^'^ and 3“^'^ nn hopping) before returning 
to the A sublattice to create a second magnon far from 
the first. Such higher energy states—ignored here but 
which we consider in the three-band model, see below— 
contribute to the spin-polaron-l-one magnon continuum 
which appears above the quasiparticle band. The rele¬ 
vance of this higher-energy feature is discussed below. 

EOM for the new propagators are generated similarly. 
We do not write them here because they are rather cum¬ 
bersome, but it is clear that the EOM for F 2 link them 


other F 2 , as well as to some of the Fi and to three- 
magnon propagators F 3 . Again we only keep those prop¬ 
agators consistent with the variational choice of having 
the 3 magnons on adjacent sites. Since 4-magnon config¬ 
urations are excluded, the EOM for F 3 link them only to 
other F 3 and to various F 2 . The resulting closed system 
of coupled linear equations is solved numerically to find 
all these propagators, including G(k, w). 

With G(k, w) known, we can find the quasiparticle dis¬ 
persion A(k) as the lowest pole of the spectral function 
A(k, w) = — i ImG(k, w). Of course, this is the quasipar¬ 
ticle in a Neel background, i.e. when the spin fluctuations 
of the AEM background are completely ignored. 


B. Quasiparticle in a background with spin 
fluctuations 

To estimate the effect of the background spin fluctu¬ 
ations (due to spin flipping of pairs of nn AFM spins, 
described by 'Hsf) we again invoke a variational princi¬ 
ple. Spin fluctuations occuring far from the ZRS should 
have no effect on its dynamics, since they are likely to 
be “undone” before the hole arrives in their neighbor¬ 
hood (they can be thought of as vacuum fluctuations). 
The spin fluctuations that influence the dynamics of the 
hole must be those that occur in its immediate vicinity 
and either remove from the quasiparticle cloud pairs of 
nn magnons generated its motion, or add to it pairs of 
magnons through such AFM fluctuations. 

For consistency, we keep the same variational config¬ 
urations here like we did at the previous step. Then, 
Eq. ([ 8 ]) acquires an additional term on the l.h.s. equal 
to: ^''"'’*'^^ 2 (k,a;, e, e'), describing pro¬ 

cesses where a pair of magnons is created through spin- 
fluctuations near the hole. Similarly, the EOMs for 
F 1 /F 2 /F 3 acquire terms proportional to F 3 /G/F 1 respec¬ 
tively, because spin fluctuations add/remove a pair of 
magnons to/from their clouds. This modified system 
of linear equations has a different solution for G(k, w), 
which accounts for the effects of the spin fluctuations that 
occur close to the hole. Comparison with the previous re¬ 
sults will allow us to gauge how important these “local” 
spin fluctuations are to the quasiparticle’s dynamics. 

Accuracy can be systematically improved by increasing 
the variational space, and implementation of such gener¬ 
alizations is straightforward. As shown next, the results 
from the variational calculation with configurations of up 
to three magnons located on adjacent sites compares well 
against available ED results and allows us to understand 
what determines the quasiparticle’s dispersion, so we do 
not need to consider a bigger variational space. 

The three-band model is treated similarly, with the 
variational space again restricted to the same configu¬ 
rations with up to three adjacent magnons. Results for 
a quasiparticle in the Neel background (no spin fluctua¬ 
tions) were published in Ref. |16l |. where the reader can 
find details about the corresponding EOMs (see also Ref. 
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(2^). Here we will focus primarily on the effect of lo¬ 
cal spin fluctuations. These are introduced as explained 
above, by adding to the EOMs terms consistent with the 
variational space and whose magnon count varies by two. 


IV. RESULTS 

A. One-band model 

We hrst present results for the one-band model. It is 
important to note upfront that it has long been known 
that both spin fluctuations and the longer range hopping 
must be included in order to obtain the correct dispersion 
for its quasiparticle [1, [Tl| . (By “correct dispersion” we 
mean one in agreement with the experimental data [^, 
0 ). Our results confirm these facts, as shown next. 

The novelty is, therefore, not in finding these results, 
but in using them to prove the validity of our variational 
method and, more importantly, to untangle the specific 
role that spin fluctuations and long-range hopping play in 
arriving at this dispersion. To the best of our knowledge, 
this had not been known prior to this work. 

The quasiparticle dispersion i?(k) is shown in Fig. [2] 
for various models: in panel (a) we set t' = t" = 0 and 
freeze spin-fluctuations {t-Jz model). In panel (b), spin 
fluctuations close to the hole are turned on as discussed; 
for simplicity we call this the t-J model, although the 
true t-J model includes all spin fluctuations. In panel 
(c), we further add the longer range hopping; for sim¬ 
plicity, we call this the model although, again, 

spin fluctuations are allowed only near the hole. Finally, 
in panel (d) we keep the longer range hopping but freeze 
the spin fluctuations; this is the t-t'-t''-Jz model. Panel 
(e) shows model dispersions explained below. 

The quality of our variational approximation is illus¬ 
trated in panels (b) and (c). Its results (thick lines) are in 
fair agreement with those of exact diagonalization (FD) 
for a 32-site cluster, which includes all spin fluctuations 
sp. Results in panel fc) agree well with those mea¬ 
sured experimentally [IflI .[llj . Our bandwidths are some¬ 
what different; some of this may be due to finite-size ef¬ 
fects, as the FD bandwidth varies with cluster size (^ . 
This also suggests that more configurations need to be 
included before full convergence is reached by our vari¬ 
ational method (these would increase the bandwidth in 
panel (b) and decrease it in panel (c), see below). This 
is supported by Ref. [1^, where full convergence for the 
t-Jz model was reached when configurations with up to 
5 magnons were included. Nevertheless, the agreement 
is sufficiently good to conclude that the essential aspects 
of the quasiparticle physics are captured by the three- 
magnon variational calculation, and to confirm that it 
suffices to include spin fluctuations only near the hole. 

These results clearly demonstrate that both spin fluc¬ 
tuations and longer-range hopping are needed to achieve 
the correct quasiparticle dispersion shown in panel (c), 
with deep, nearly isotropic minima at (^,^). Absence 



FIG. 2. (color online) Quasiparticle energy E(k) along several 
cuts in the Brillouin zone for various one-band models. In all 
cases J/t — 0.3 while t' = t" — 0 in (a),(b), and t'/t = 
—Q.3,t"/t = 0.2 in (c),(d). Lines show the results of the 
variational calculation with the spin fluctuations frozen in (a) 
and (d), or allowed only near the hole in (b) and (c). Symbols 
in (b) and (c) are the corresponding ED results for a 32- 
site cluster [5 0. (e) Dispersion iJgp(k) of Eq. (fTTIl for 

Eo = 0 and t 2 = —l,t 3 = 0 (black); t 2 = l,t 3 = 0.5 (red); 
t 2 = -l;t3 = 2/3 (green). 


of longer-range hopping leads to a rather flat dispersion 
along (ffi'^) — (tt, 0 ), see panel (b); this fact has long been 
known [9|. In panel (d), we show that if the longer range 
hopping is included but spin fluctuations are absent, the 
dispersion is rather flat along the (0,0 ) — (tt, tt) direction. 
We are not aware of previous studies of this case. 

Although F(k) looks very different in the four cases, it 
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FIG. 3. (color online) (a) Shortest Trugman loop that generates a t 2 ,TL contribntion. A mnch smaller is also generated, 
see Ref. [^; Generation of (b) t 2 ,sf and (c) terms due to spin-fluctuations; (d) Process that renormalizes t". The square 
shows the location of the hole, while the circles show magnons (wrongly oriented spins). The remaining spins are in their Neel 
order orientation and are not shown explicitly. The short thick arrows indicate the next step in the process, while the thin 
arrows in the final sketch show the effective quasiparticle hopping generated by those processes. 


turns out that all can be understood in a simple, unified 
picture. The key insight is that the quasiparticle lives on 
one magnetic sublattice, because of spin conservation. As 
a result, its generic dispersion must be of the form: 

EqpCk) = Eq + 4t2 cos kx cos ky -I- 2t3(cos 2kx + cos 2ky) 

( 11 ) 

i.e. like the bare hole dispersion e(k) of Eq. Q, but with 
renormalized 2'^'^ and 3'^'^ nn hoppings t' —>■ t 2 ; t” —>■ ^3 ■ 
There cannot be any effective nn hopping of the quasi¬ 
particle because this would move it to the other sublat¬ 
tice; this cannot happen without changing the magnetic 
background, so ti = 0. Longer range hoppings that keep 
the quasiparticle on the same sublattice may also be gen¬ 
erated dynamically, but their magnitude is expected to 
be small compared to ^ 2 ,^ 3 , hence Eq. m- Thus, un¬ 
derstanding the shape of the quasiparticle dispersion re¬ 
quires understanding the values of t 2 and 

We begin the analysis with the t-J^ model. Its quasi¬ 
particle is extremely heavy, as shown in panel (a). Note 
that the vertical scale is an order of magnitude smaller 
than for the other panels. The reason is that every time 
the hole hops, it moves to the other magnetic sublattice 
and it must either create or remove a magnon, to con¬ 
serve the total spin. As the hole moves away from its 
original location, it leaves behind a string of magnons 
whose energy increases roughly linearly with its length. 
This could be expected to result in confinement (infinite 
effective mass), but in fact the quasiparticle acquires a 
finite dispersion by executing Trugman loops (TL) [^ . 
the shortest of which is sketched in Fig. EJa). By go¬ 
ing nearly twice along a closed loop, creating a string 
of magnons during the first round and removing them 
during the second round, the hole ends up at a new lo¬ 
cation on the same magnetic sublattice. Only 2'^'^ and 
nn hopping terms can be generated through TL ir¬ 
respective of their length, and 1 ^ 3 , tl| "C IO.tlI -C T if 
t/J ~ 3 [ 2 ^. Indeed, setting 0 < 0,t3 —>■ 0 in ifqp(k) of 
Eq. (fTTl) leads to the black curve in Fig. [2Ke) [ 2 ^, which 
has the same shape as that of panel (a) (the bandwidth 
is proportional to |t 2 ,TL|)- This dispersion is wrong not 
just quantitatively but also qualitatively, with (f, f) as 
a saddle point instead of the ground state. Clearly, ig¬ 
noring both longer range hopping and spin fluctuations 
changes completely the dynamics of the quasiparticle. 


When the spin fluctuations are turned on in the t-J 
model, they act on a time scale Tgf ^ 1 /J much faster 
than the slow dynamics due to TL, ttl ~ 1 / 1 ^ 2 , tl|- The 
main contributions to t 2 and now come from processes 
like those sketched in Fig. [3Kb) and (c), where spin fluc¬ 
tuations remove pairs of magnons created by nn hop¬ 
ping of the hole, leading to t 2 .sf > |i 2 ,TL|, ts.sf > I^s.tlI- 
Moreover, we expect t 2 ,sf = 2 t 3 _sf because these effective 
hoppings are generated by similar processes but there are 
twice as many leading to 2 "^^ compared to 3'''^ nn hopping, 
as the hole can move on either side of a plaquette. 

Indeed, the t-J dispersion of panel (b) has a shape sim¬ 
ilar to that of ifnn(k) with t 2 = 2 t 3 , shown as a red curve 
in Fig. j^Je) [s^. Because Eqp{k,TT — k) = Eq — 2 t2 + 
( 4 t 3 — 212 ) cos 2 fc, it has a perfectly flat dispersion along 
(0,7r) — (tt, 0) for t 2 = 2t^. The dispersion in panel (b) 
is not perfectly flat along this cut, so in reality ^2 ~ 2 t 3 . 
The small correction from the factor of 2 is likely due 
to higher order processes, as well as contributions from 
TL (which remain active). Ignoring it, we find the corre¬ 
sponding bandwidth Eqp[Q, 0 ) — Eqp{^, ^) = 4^2 + 8^3 = 
8 ^ 2 , suggesting that the effective hoppings generated with 
spin fluctuations are of the order t 2 ,sf ~ 2t3.sf ~ J/E 

Next, we consider what happens if instead of (local) 
spin fluctuations, we turn on longer-range hopping. Un¬ 
like in the t-Jz model, the quasiparticle of the t-t'-t”- 
Jz model should be light because the longer range hop¬ 
pings allow the hole to move freely on its magnetic 
sublattice. It can therefore efficiently remove magnons 
created through its nn hopping, without having to com¬ 
plete the time-consuming Trugman loops. The presence 
of the magnon cloud renormalizes these bare hoppings to 
smaller values, as is typical for polaron physics. Figure 
|3Kd) shows one such process that renormalizes t” t"*. 
Similar processes (not shown) renormalize t' —> t'* so 
both hopping integrals should be renormalized by com¬ 
parable factors. As a result, we expect a dispersion like 
Eqp{k) but now with ^ 2/^3 = t'*It"* k, t' jt", if we ig¬ 
nore the small TL contributions. This indeed agrees with 
the result in panel (d), as shown by its comparison with 
the green curve in Fig. H e) where Fqp(k) is plotted for 
^ 2/^3 = t'jt" = —1.5. For t 2 = —2^3, Fqp(k) would be 
perfectly fiat along (0,0) — (tt, tt). Thus, the change in the 
relative sign explains why now the dispersion is nearly fiat 
along (0,0) — (tt, tt) and maximal along (0, tt) — (tt, O), in 








contrast to the previous case. However, while t 2 ,st/i 3 ,sf ~ 
2 is always expected for the t-J model so its dispersion 
must have a shape like in Fig. [SJb), in the 
model the ratio t'*/t"* mirrors the ratio t'/t". If this 
had a very different value than Ri —2, the quasiparticle 
dispersion would change accordingly. 

These results allow us to understand the dispersion of 
the quasiparticle. This must have contributions 

from both the spin fluctuations and the renormalized 
longer range hoppings, plus much smaller TL terms, be¬ 
cause the processes giving rise to them are now all active. 
Indeed, the curve in panel (c) of Fig. [2]is roughly equal to 
the sum of those in panels (b) and (d). The isotropic min¬ 
imum at (f, f) is thus an accident, since the dispersion 
along (0,0) — (tt, tt) is controlled by spin fluctuations, and 
that along (0, tt) — (tt, 0) is due to the renormalized longer 
range hoppings. More precisely, because t 2 ,sf ~ 
the contributions coming from spin fluctuation interfere 
destructively for momenta along (0, tt) — (tt, 0) so disper¬ 
sion here is controlled by the renormalized t'* k, — 1.5t"*, 
and viceversa. If t 2 ,sf ~ (which happens to hold be¬ 
cause J ^ |t'|), the sum gives nearly isotropic dispersion 
near (f,f). If we change parameters significantly, the 
dispersion becomes anisotropic (not shown). 

Before moving on to contrast this behavior with that of 
the quasiparticle of the three-band model, we briefly dis¬ 
cuss the effect of the three-site term of Eq. ©■ The vari¬ 
ational results for the four models are shown in Fig. S) 
Where direct comparisons can be made, they are again in 
good quantitative agreement with other work where this 
term has been included, such as in Ref. [3l|. Its inclusion 
has a qualitative effect only for the t-J^ model, where the 
shape of the dispersion is changed in its presence. This 
is not very surprising because, as discussed, the Trug- 
man loops which control behavior in that case are very 
slow processes, and their effect can easily be undone by 
terms that allow the hole to move more effectively. The 
three-site term is such a term and its presence increases 
the bandwidth not just for the t-Jz model, but for all 
cases. For the other three models, however, the inclusion 
of this term changes the dispersion only quantitatively: 
the bandwidth is increased but the overall shape is not 
affected much. The biggest change is along (0, 0) — (tt, tt), 
as expected because the three-site term generates effec¬ 
tive 2nd and 3rd nn hoppings with the same sign and a 
2/1 ratio, i.e. similar to t 2 ,s/ and As a result, its 

presence mimics (and boosts) the effect of the local spin 
fluctuations. 

It is interesting to note that if we allow this term to 
be large enough, we could obtain a dispersion with the 
correct shape even in the absence of spin fluctuations. 
However, the scale of this term is set by J, it is not a free 
parameter. As a result, we conclude that with its proper 
J energy scale, this terms does not change qualitatively 
the behavior of the quasiparticle of the one-band model 
(apart from the t-Jz case), although its inclusion may, in 
principle, allow for better fits of the experimental data. 



_I_I_I_^_I_I_I_ 

( 71 , 0 ) ( 0 , 0 ) ( 71 , 71 ) ( 0 , 71 ) ( 71 , 0 ) 


FIG. 4. (color online) Quasiparticle dispersion when the 
three-site term of Eq. © is included in the one-band Hamil¬ 
tonians (dashed lines). For comparison, the dispersions with¬ 
out this term are also show (full lines from Fig. [2| 


B. Three-band model 

Results for the simplified three-band model with the 
spin-fluctuations frozen off were discussed in Ref. M- 
To keep this work self-contained, we show in Fig. [5] 
the most relevant data for the issue of interest, namely 
the quasiparticle dispersion Efk) obtained in a varia¬ 
tional calculation with the maximum number of magnons 
n-m = 0 — 3. These results already suffice to illustrate the 
qualitative difference between the quasiparticle dynamics 
in the one-band and the three-band models. 

The Um = 0 curve plots the dispersion if no magnons 
are allowed, i.e. not only the spin fluctuations of the 
AFM background but also spin-flip processes due to Jpd 
and Tswap are turned off. It is important to emphasize 
that the resulting dispersion does contain a very impor¬ 
tant contribution from the terms in Tswap describing hop¬ 
ping of the hole past Cu spins with the same spin projec¬ 
tion, so that the spin-swap leaves the spins unchanged. 

In fact, it is the interference between these terms with 
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FIG. 5. (color online) E(k) along several cuts in the Bril- 
louin zone for the three-band model. The results are for the 
variational calculation with the spin-fluctuations turned off 
and conhgurations with up to Um magnons allowed. The 
“restricted” calculations labelled n 2 ,r,n 3 ,r imposed the ad¬ 
ditional constraint that the magnons are on adjacent sites. 
While the bandwidth is strongly renormalized with increasing 
Um, the nearly isotropic dispersion around the ground-state 
at (f, f) is a consistent feature. See text for more details. 


those of Tpp that leads to this interesting bare disper¬ 
sion, which already has deep, nearly-isotropic minima 
near (f, f). If we set tsw = 0, the bare dispersion due 
to only Tpp has the ground-state at (7r,7r), whereas if 
Tpp = 0, the dispersion due to the allowed terms in Tg^ap 
is perfectly flat because the hole is then trapped near a 
like Cu spin. However, as long as tpp ~ tsw, the isotropic 
minimum emerges at (^,^). In this context, it is use¬ 
ful to note that in many numerical studies of the Emery 
model, tpp was set to zero simply for convenience. Our 
results suggest that this choice changes the quasiparticle 
dynamics qualitatively, and is therefore unjustified. 

The bare hole dispersion in the three-band model thus 
already mimics this key aspect of the correct quasipar¬ 
ticle dispersion, unlike in the one-band model. Allowing 
the quasiparticle cloud to emerge by allowing the hole 
to create and absorb magnons in its vicinity, through 
spin-flip processes controlled by Jpd and T’g.^yap) further 
renormalizes the bandwidth (a typical polaronic effect) 
without affecting the existence of the isotropic dispersion 
near (^, ^). This magnon cloud is very important, how¬ 
ever, to stabilize the low-energy quasiparticle, as demon¬ 
strated by the significant lowering of the total energy. In 
particular, at least one magnon must be present in order 
for a ZRS-like object to be able to form, and indeed the 
TLm = 1 curve is pushed down by ^ \0Jdd compared to 
the bare dispersion. We further analyze the relevance of 
the ZRS solution below. 

The small difference between the rim = 2 and rim = 2, r 
results proves that magnons indeed sit on adjacent sites 
in the cloud. (The latter solution imposes this constraint 


explicitly, whereas the former allows the magnons to be 
at any distance from each other. In both cases, the hole 
can be arbitrarily far from the magnons, although, as ex¬ 
pected, configurations where the hole is close to the last 
emitted magnon have the highest weight in the quasi¬ 
particle eigenstates.) At higher energies, however, these 
two solutions are qualitatively different. The former con¬ 
tains the expected quasiparticle-l-one-magnon continuum 
starting at Ei ^s + 2Jdd, where Ei ^s is the ground-state 
energy of the quasiparticle with rim = Ij and 2 Jdd is the 
energy cost to create a magnon far from it. Their sum 
is the energy above which higher-energy (excited) states 
must appear in the spectrum, describing the quasiparti¬ 
cle plus one magnon not bound to its cloud. The presence 
of this continuum guarantees that in the fully converged 
limit, the quasiparticle bandwidth cannot be wider than 
2 Jdd, since the quasiparticle band is always “flattened 
out” below this continuum (another typical polaronic be¬ 
havior). For both n 2 ,r and calculations, the quasi¬ 
particle is already heavy enough that its dispersion fits 
below the corresponding continuum. This is why enlarg¬ 
ing the variational space with configurations needed to 
describe this feature, with at least one magnon located 
far from the cloud, does not affect the quasiparticle dis¬ 
persion much (see Ref. [l^ for more discussion). 

The bandwidth of the rim = 3, r dispersion is in de¬ 
cent agreement with numerical results for this model, as 
discussed next, suggesting that this variational calcula¬ 
tion is close to fully converged. The fact that the cloud is 
rather small should not be a surprise. The variational ap¬ 
proach explicitly imposes the constraint that there is at 
most one magnon at a site. As magnons sit on adjacent 
sites when bound in the quasiparticle cloud, they prefer 
to occupy a compact area to minimize their exchange en¬ 
ergy cost, thus creating a domain in the other Neel state 
(down-up instead of up-down). The hole prefers to sit 
on the edge of this domain, because being inside it is 
equally disadvantageous to being outside, i.e. far from 
magnons. However, since on the boundary the hole can 
interact with only one magnon at one time, a large and 
costly domain is unlikely. 

We can now contrast the dynamics of the quasiparticle 
in the three-band model if the spin fluctuations are frozen 
out with the corresponding one-band model, namely the 
t-t'-t''-Jz case. Both have a quasiparticle with a small, 
few-magnon cloud, and a bandwidth k. 2J = 2Jdd- The 
key difference is that the three-band model already shows 
a dispersion with the correct shape, whereas for the one- 
band model the dispersion is much too flat along (0,0) — 
(7r,7r). This difference is traced back to the fact that in 
the one-band model, the bare hole dispersion also suffers 
from this same problem if t'/t" ^ —1.5, unlike that of 
the three-band model. As a result, spin fluctuations are 
necessary to find the correct dispersion in the one-band 
model, as already shown, but their role in the three-band 
model should be rather limited. 

To confirm this conjecture, we consider the effect of lo¬ 
cal spin fluctuations on the dispersion of the three-band 
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FIG. 6. (color online) i?(k) along several cuts in the Brillouin 
zone for the three-band model in the restricted variational ap¬ 
proximations with (a) Um = 2 and (b) rim = 3. Circles show 
ED results for St = | from Ref. for a 32 Cu + 64 O clus¬ 
ter, shifted to have the same ground-state energy. Full lines 
show the results of Fig. O without spin fluctuations. Orange 
lines with square symbols are the results if spin fluctuations 
occur near the hole. The dashed green line in panel (b) is 
the dispersion when spin fluctuations are allowed to locally 
create/remove a pair of magnons only if no other magnons 
are present/remain in the system. See text for more details. 

model quasiparticle. In Fig. [6l results for the restricted 
variational approach with rim = 2,r and rim = 3,r [i.e. 
up to two or up to three magnons on a^acent sites) are 
compared to the ED results of Ref. [^, shown by the 
black full circles. Full lines (red and blue, respectively) 
show the results of Fig. [3 without spin fluctuations. 
Orange lines with squares show the dispersion with spin 
fluctuations turned on near the hole. The dashed green 
line in panel (b) shows an intermediate result when spin 
fluctuations are allowed to create a pair of magnons only 
if there is no magnon in the system, and to remove a pair 
if only two magnons are present (the orange line also in¬ 
cludes contributions from processes where spin fluctua¬ 
tions add a pair of magnons when a magnon is already 
present, and its reversed process). 

The effect of spin fluctuations is similar to that found 
in the one-band models, as expected because the AFM 
background is modeled identically. They again have very 
little effect on the (tt, 0) — (0, tt) dispersion; beside a small 
shift to lower energies, this bandwidth is only slightly in¬ 
creased, bringing it into better agreement with the ED 
values for rim = 3. Like for one-band models, spin 
fluctuations lead to a more significant increase of the 
(0,0) — (7r,7r) dispersion. For rim = 3, it changes from 
being too narrow without spin fluctuations, to too wide 
in their presence. (The rim = 2 overestimate of the band¬ 
width is expected, see discussion in Ref. 0 )- 

The increased energy near (0, 0) = (tt, tt) may seem 
problematic but one must remember that in reality, F(k) 
is flattened below a continuum that appears at 
above the ground state. The continuum is absent in 


this restricted calculation because configurations with a 
magnon far from the cloud, which give rise to it, are not 
included. This explains why the overestimated band¬ 
width is possible. In the presence of the continuum, 
states that overlap with it hybridize with it and a dis¬ 
crete state (the quasiparticle) is pushed below its edge. 
This will lower the value at (0,0) and lead to good quan¬ 
titative agreement everywhere with the ED results. (The 
keen reader may note that the ED bandwidth is also 
slightly wider than but one must remember both 

the finite size effects of cluster ED, and the existence of 
a S't = I polaron in its spectrum Since the total 

spin Sy is not a good quantum number in our variational 
approximation, its quasiparticle probably overlaps some¬ 
what with both the St = \ and St = \ spin-polarons, 
so it is not clear which ED states to compare against). 

Our results show that spin fluctuations have a similar 
effect in both models. However, while they are essential 
for restoring the proper shape of the dispersion in one- 
band models, they are much less relevant for the three- 
band model. This is a direct consequence of the different 
shape of the bare bands, as discussed, but also to having 
J ~ |t'| while Jdd ~ tpp/A. In the three-band model, the 
quasiparticle creates and absorbs magnons while moving 
freely on the O sublattice, on a timescale that is faster 
than that over which spin fluctuations act, and so their 
effect is limited. In contrast, in the one-band model, the 
timescale for free propagation of the hole on the same 
magnetic sublattice (controlled by t',t'') is comparable 
with the spin fluctuations’ timescale, and therefore the 
effect of spin fluctuations is much more significant. They 
are especially important along (0,0) — (7r,7r), where the 
bare dispersion of one-band models is nearly flat. 


V. DISCUSSION AND SUMMARY 

In this work, we used a variational method to study 
and compare the quasiparticle of t-J and one- 

band models, to that of a (simplified) three-band model 
that is the intermediary step between the full three-band 
Emery model and the one-band models. 

Our variational method generates the BBGKY hierar¬ 
chy of equations of motions for a propagator of interest 
(here, the retarded one-hole propagator), but simplified 
by setting to zero the generalized propagators related 
to projections on states that are not within the varia¬ 
tional space. Its physical motivation is very simple: if 
the variational space is properly chosen, i.e. if it con¬ 
tains the configurations with the highest weight contri¬ 
butions to the quasiparticle eigenstates, then the ignored 
propagators are indeed small because their residue at the 
Lj = £'(k) pole is proportional to their weight (Lehmann 
representation). Setting them to zero should thus be an 
accurate approximation. Numerically, the motivation is 
also clear: because the resulting simplified hierarchy of 
coupled equations can be solved efficiently, we can quite 
easily study a quasiparticle (or a few [s^) on an infinite 
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plane, thus avoiding finite-size effects and getting full in¬ 
formation about k dependence, not just at a few values. 
Moreover, by enlarging the variational space and by turn¬ 
ing off various terms in the Hamiltonian, both of which 
lead to changes in the EOM and thus the resulting prop¬ 
agators, one can infer whether the calculation is close to 
convergence and isolate and understand the effect of var¬ 
ious terms, respectively. The ability to efficiently make 
such comparisons is essential because it allows us to gain 
intuition about the resulting physics. 

Our results show that even though for reasonable val¬ 
ues of the parameters, the quasiparticle dispersion £l(k) 
has similar shapes in both models, the underlying quasi¬ 
particle dynamics is very different. In the three-band 
model, the bare dispersion of the hole on the O sub¬ 
lattice, due to tpp and spin-swap hopping tsw past Cu 
with parallel spins, already has a deep isotropic minimum 
near (f, f), unlike the bare £(k) of the one-band models. 
When renormalized due to the magnon cloud, it produces 
a quasiparticle dispersion with the correct shape in the 
whole Brillouin zone, even in the absence of spin fluctua¬ 
tions. In contrast, for the one-band models the inclusion 
of spin fluctuations is necessary for the correct dispersion 
to emerge. This shows that the quasiparticle dynamics 
is controlled by different physics in the two models, and 
this is likely to play a role at finite concentrations as well. 

Our results thus raise strong doubts on whether the 
one-band model truly describes the same physics 

like the three-band model. We can think of three possible 
explanations to explain these differences: 

(i) The model is the correct one-band model, 

but its true parameters have values quite different from 
the ones used here. Indeed, if the bare e(k) dispersion 
had isotropic minima at (§, ^), and if its renormalized 
bandwidth would be of the order of 2J, then spin fluc¬ 
tuations could not change it much, similar to what is 
observed in the three-band model. 

This explanation can be ruled out. An isotropic bare 
dispersion e(k, k) « e(fc, n — k) requires that \t''/t'\ » 1, 
which is physically unreasonable. 

(ii) The model has a different quasiparticle 

because its underlying assumption, i.e. the existence of 
the low-energy ZRS, is wrong. This would mean that 
not only this specific one-band model but any other one 
obtained through such a projection would be invalid. 

We can test this hypothesis by calculating the overlap 
between the three-band quasiparticle and a ZRS Bloch 
state. The latter is defined in the only possible way that 
is consistent with Neel order: 


|ZRS.k) = — ^ 
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where 

A 




Pl+^,a + Pl+^,a - - pI-^. 


is the linear combination of the p orbitals neighbor to the 
Cu located at i, that has the overall x'^ — symmetry 



FIG. 7. (color online) (a) Overlap pzRS between the ZRS 
Bloch state of Eq. m and the quasiparticle eigenstate, 
as obtained in the restricted variational calculations with 
rim = 1,2,3, without (full lines) and with (full symbols) lo¬ 
cal spin fluctuations included. The empty squares and the 
dashed line show the weight of the ZR triplet, (b) Overlap 
PzRS normalized with respect to the probability to have no 
magnon (po) or to have one magnon near the hole (pi), in the 
quasiparticle eigenstate. See text for more details. 


(our choice for the signs of the lobs is shown in Fig. [1}. 
Note that with this definition |(ZRS,k|ZRS, k)p = I for 
any k, so there are no normalization problems Q. 

We define pzRS = |( 9 T,k|ZRS,k)p as the overlap be¬ 
tween the quasiparticle eigenstate of momentum k and 
this ZRS Bloch state. Its value can be calculated from the 
appropriate residues of the zero- and one-magnon prop¬ 
agators at a; = F(k), and is shown in Fig. [TKa). We do 
not plot the rim = 0 results because a singlet cannot form 
if the Cu spins cannot flip. (For rim = 0, there is over¬ 
lap with the spin-up hole component of |ZRS, k), and we 
find that pzRS varies from 0 at (0, 0) and (0,7r) to 0.5 at 
(^, ^), but the same answer is found for a triplet. Inter¬ 
estingly, this proves that the bare hole dispersion already 
has eigenstates with the symmetry near (f, f)). 

For rim = 1 "we find pzRS 0-9 in the entire Brillouin 
zone. Clearly, in this very small variational space, locking 
into a ZRS is the best way for the doped hole to lower 
its energy. However, the value of pzRS decreases fairly 
significantly for rim = 2,7- and rim = 3, r. First, note 
that turning the spin fluctuations on or off has almost 
no effect on pzRS ■ This is consistent with our conclusion 
that local spin fluctuations do not influence the nature 
of the quasiparticle in the three-band model: clearly, its 
wavefunction is not changed in their presence. 

The decrease of pzRS with increasing rim could be 
due either to increased contributions to the eigenstate 
from many-magnon configurations (which have no over¬ 
lap with |ZRS,k)), and/or from competing states such 
as a ZR triplet, and/or singlets or triplets with the hole 
occupying a linear combination of O orbitals with s,px 
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or py instead of symmetry. The latter possibility 

can also be ruled out because overlaps with those Bloch 
states are found to be small. The largest such contri¬ 
bution is from the ZR triplet state, shown in Fig. [Tl^a) 
by the dashed line and open squares for rim = 3, r with¬ 
out and with local spin fluctuations, respectively. This 
overlap is much smaller than with the ZRS singlet. 

Another way to confirm this is displayed in Fig. [71[b), 
where we compare pzRS to pQ+pi, where po is the prob¬ 
ability to find the hole without any magnons, and pi is 
the probability to find one magnon adjacent to the hole, 
in the quasiparticle eigenstate. Note that Po + Pi < 1 
even for the one-magnon variational approximation be¬ 
cause the hole can also be located away from the magnon. 
As rim increases, po -I- pi decreases even more as config¬ 
urations with two or more magnons now also contribute 
to the normalization. These configurations with two or 
more magnons, and those with one magnon not adjacent 
to the hole, have no overlap with |ZRS,k), explaining 
the decrease in the magnitude of pzRS- However, the 
ratio pzRs/(Po +Pi) > 0-9 in the whole Brillouin zone, 
confirming that this part of the wavefunction has a pre¬ 
dominant ZRS-like nature. This is certainly the case near 
the (f, f) point, where the overlap is converged to 1. In¬ 
terestingly, at the antinodal points this ratio decreases 
with increasing rim, and here the overlap with the ZR 
triplet is largest, see Fig. Ela), suggesting that a ZRS 
description is less accurate in this region. 

Thus, the zero- and one- (adjacent) magnon parts of 
the wavefunction have significant overlap with the ZRS 
Bloch state. However, pzRS ~ 0.5 is a rather small 
value, and it is not clear whether the dressing with 
more magnons is consistent with this ZRS picture or not. 
It is possible that the two- and three-magnon compo¬ 
nents of the wavefunction have significant overlap with a 
ZRS-bone magnon and ZRS-l-two magnon configurations, 
but they could also have quite different nature. It is not 
clear to us how to verify which is the actual situation. 

If these two- and three-magnon components have sig¬ 
nificant non-ZRS character, however that is defined in 
this case, then clearly the difference observed in the re¬ 
sults from one- and three-band models would be likely 
due to this non-ZRS nature. 

If, on the other hand, one takes these results to support 
the idea that a low-energy projection onto ZRS states is 
valid, then this is not the origin of the discrepancy in the 
quasiparticle behavior. In this case, it must follow that: 

(hi) The is not the correct one-band model be¬ 

cause there are additional important terms generated by 
the projection onto the ZRS states, like those discussed 
in Ref. or the three-site terms, which it neglects. 

If (iii) is indeed the explanation for the different be¬ 
havior of the quasiparticles of the one- and three-band 
models, then in our opinion this implies that the strategy 
of using one-band models to study cuprates is unlikely to 
succeed. The main reason for this strategy, as mentioned, 
is to make the Hilbert space as small as possible for com¬ 
putational convenience. This, however, is only useful if 


the Hamiltonian is also fairly simple. 

In principle one could test additional terms that could 
be included in one-band models by using methods like 
ours, to figure out which insure that the resulting be¬ 
havior mirrors that of the three-band model. Even if 
this enterprise was successful and the “fix” was relatively 
simple, i.e. only a few additional terms and correspond¬ 
ing parameters are necessary, it is important to empha¬ 
size that this improved one-band model would still not 
describe correctly cuprates at finite doping. Additional 
terms must be included to correctly account for the ef¬ 
fective interactions between quasiparticles in one-band 
models, as we demonstrated in Ref. [s^. 

From a technical point of view, their origin is simple 
to understand. Even for the simplified three-band model, 
the presence of additional holes leads to additional terms 
in the Hamiltonian [^ . because the intermediary states 
are different and this affects the projection onto states 
with no-double occupancy on Cu. This is going to be¬ 
come even more of an issue if a subsequent projection 
onto ZRS-like states has to be performed, and may well 
result in an unmanageably complex Hamiltonian. 

This is why we believe that the (simplified) three-band 
model is a safer option to pursue. Its computational com¬ 
plexity is not that much worse than for one-band models, 
whereas the Hamiltonian is certainly simpler. In fact, our 
demonstration here that there is no need to accurately 
capture the spin fluctuations of the AEM background in 
order to gain a reasonable understanding of the quasi¬ 
particle behavior, makes its study significantly simpler. 
In particular, it allowed us to study one hole on an infi¬ 
nite layer very simply and efficiently. Generalizations to 
few holes [S^ and to finite concentrations could also turn 
out to be easier to carry out than the effort of finding the 
correct form for a one-band Hamiltonian. 

Of course, there is no guarantee that the (simplified) 
three-band model captures all physics needed to explain 
cuprates, either. It is possible that important aspects of 
the Emery model were lost through the projection onto 
spin degrees of freedom at the Cu sites (this, however, 
would affect the one-band models just as much). Even 
the Emery model itself may not be general enough; for 
instance, a generalization to a 5-band model including 
non-ligand O 2p orbitals might be needed, as suggested 
recently in Ref. . We note that such a generalization 
can be easily handled by our method (provided that one 
can still project onto spin degrees of freedom at the Cu 
sites), as showed in Ref. [l^ where we found that these 
states do not change the quasiparticle dispersion much, 
although they do have an effect on its wavefunction. 

While careful investigation of such scenarios is left as 
future work, one clear lesson from this study is that ob¬ 
taining the correct dispersion for the quasiparticle of an 
effective model is not sufficient to validate that model. 
The dispersion can have the correct shape for the wrong 
reasons, as we showed to be the case for the t-t'-t''-J 
model, where it is due to the interplay between the ef¬ 
fects of the longer-range hopping and spin-fluctuations. 
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The same dispersion is obtained for the simplified three- lent despite the similar dispersion of their quasiparticles, 
band model, however in this case the spin-fluctuations 
play essentially no role, so the underlying physics is very 

different. This difference is very likely to manifest itself in ACKNOWLEDGMENTS 
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